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Governing equation? of motion are derived for a bulk-flow model of 
the leakage path between an impeller shroud and a pump housing. The 
governing equations consist of a path-momentum, a circumferential - 
momentum, and a continuity equation. The fluid annulus between the 
impeller shroud and pump housing is assumed to be circumferentially 
symmetric when the impeller is centered; i.e., the clearance can vary 
along the pump axis but does not vary in the circumferential direction. 

A perturbation expansion of the governing equations in the 
eccentricity ratio yields a set of zeroth and first-order governing 
equations. The zeroth-order equations define the leakage rate and the 
circumferential and path velocity distributions and pressure 
distributions for a centered impeller position. The first-order 
equations define the perturbations in the velocity and pressure 
distributions due to either a radial-displacement perturbation or a 
tilt perturbation of the impeller. Integration of the perturbed 
pressure and shear-stress distribution acting on the rotor yields the 
reaction forces and moments acting on the impeller face. 

Introduction 

Figure 1 illustrates an impeller stage of a multi-stage 
centrifugal pump. Leakage along the front side of the impeller, from 
impeller discharge to inlet, is restricted by a wear-ring seal, while 
leakage along the back side is restricted by either an interstage seal 
or a balance-discharge seal. Lomakin [1] originally recognized the 
major influence that seals have on the rotordynamic response of 
centrifugal pumps. Recent analysis and test results are provided by 
references [2-4], 

More recently, various investigators have considered the forces 
developed by flow through the impeller and its interaction with either 
a volute or a vaned diffuser. Cal Tech researchers [5] have presented 
measured force coefficients for an impeller precessing in a volute. 
Ohashi and Shoj i [6] also provide measured force coefficients for an 
impeller whirling in vaneless and vaned diffusers. More recently, 
Bolleter et al. [7] from Sulzer Brothers, Ltd. have also presented test 
results for an impeller in a vaned diffuser. The Cal Tech and Sulzer 
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test programs use conventional impellers in their test programs, but 
use a radial face seal to minimize or eliminate the forces which would 
normally be developed by the wear-ring seals. The face seals employed 
by these investigators sharply reduce the leakage flow down the front 
face of the impeller. The impeller tested by Ohashl et al. has flat, 
parallel, front and back shrouds. Test results from all these 
investigators suggest that impellers create relatively benign forces 
from a rotordynamic viewpoint. Ohashi concludes that impeller forces 
would damp forward whirling motion in most operating conditions. The 
Cal Tech and Sulzer researchers conclude that impellers would provide 
positive damping unless a rotor's running speed exceeded its critical 
speed by an approximate factor of 2.5. Hence impellers are predicted 
to yield smaller destabilizing forces than bearings or long seals which 
are predicted to yield destabilizing forces when the running speed 
exceeds twice the critical speed. 

There is some contrary evidence from operating pump experience 
with respect to impeller forces. Specifically, Massey [8] cites 
experience with an 11 -stage pump which was unstable and whirled at 80$ 
of running speed; i.e., the pump was unstable when its running speed 
was only 1.25 times the critical speed. The HPOTP (High Pressure 
Oxygen Turbopump) of the SSME also whirled at 80$ of running speed [9]. 
This evidence suggests that some unaccounted-for destabilizing force is 
present in pumps. The present analysis is aimed at investigating the 
forces and moments developed by impeller-shroud forces. 

A bulk-flow analysis is employed similar to that of Childs and Kim 
[3] . However the analysis is extended to account for the changing 
geometry of impeller surface. In addition, the shear stress 
contribution to the reaction forces is accounted for, and the reaction- 
moment coefficients are calculated using the approach of reference 
[10]. Specifically, the following general model is used to define the 
reaction forces and moments which arise for small motion of a pump 
impeller about its centered position. 
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where (Fy, Fy), (X, Y) define the components of the impeller reaction 
forces and relative displacements, and (My, My), (ay, ay) define the 
components of the reaction moments and rotation (small-angles) vectors. 

Geometry and Kinematics 

Figure 1 illustrates the annular leakage paths along the front and 
back sides of a typical shrouded impeller of a multistage centrifugal 
pump. The present discussion concentrates on the flow and pressure 
fields within the forward annulus; however, the analysis also applies 
to the rear annulus. As illustrated in figure 2, the outer surface of 
the impeller is a surface of revolution formed by rotating the curve 
R - R(Z) about the Z axis. A point on the surface may be located by 
the coordinates Z, R(Z), 0. The length along the curve R(Z) from the 

initial point R^ , Z ± to an arbitrary point R, Z is denoted by S and 
defined by 

Z _ R 

S - / / 1+ dR 2 du « / / 1 + dZ 2 du (2) 

Zy dZ Rj_ dR 

In the equations which follow, the path coordinate S and angular 
coordinate 0 are used as independent spatial variables. The co- 
ordinates Z,R defining the impeller surface are be expressed as 
parametric functions of S, i.e., Z(S), R(S). 

Trigonometric functions of the angle Y , illustrated in figure 3, 
are defined as follows 


dR 


tanY 

= - — 




dZ 


(3) 


dZ 

dR 


cosY 

«= — , sinY « - 

— 



dS 

dS 


The clearance between 

the impeller 

and 

the housing is denoted as 


H(S ,0 , t) , with the time dependency introduced by impeller motion. In 
the centered position, the clearance function depends only on S and is 
denoted by H 0 (S). Displacement of the impeller in the X and Y 
directions obviously causes a change in the clearance function. The 
clearance function is also changed by pitching or yawing of the 
impeller as defined by rotations ay, ay about the X and Y axes, 
respectively. For small displacements and rotations of the impeller 
the clearance function can be stated 

H(S,0,t) - H 0 (S) - [(X + ayZ) cosY - ayRsinY] cos0 

(4) 

“ [(Y - ayZ) cosY + ayRsinY] sin0 

Observe in this equation that H 0 , R, Z, cosY and sinY are solely 
functions of S, while X, Y, ay, ay are functions only of t. 
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Governing Equations 

Returning to figure 2, the path coordinate S and circumferential 
coordinate R6 are used to locate a fluid differential element of 
thickness H(S,6,t) illustrated in figure 4. From the geometry of 
figure 4 , the continuity equation can be stated 

3H 3 13 H 3R 

— + — (U S H) + (U 0 H) + (-) — U s - 0 

3t 3S R 36 R 3S 

where U s and Uq are the path and circumf erential bulk-velocity 
components, respectively. 

Figure 5 illustrates the pressure and shear-stress components 
acting on the differential fluid element. The first subscripts (s,0) 
in the shear-stress definitions (t ss , T sr ), (xe ? , Tg r ) denote path and 
6 directions, respectively; the second subscripts (s, r) denote stator 
and rotor surfaces, respectively. The path momentum equation can be 
stated 


3P U e 2 dR 

-H — - - pH + t S5 + t sp 

3S R dS 


3U 


+ PH (— * 
3t 


s 3U S Uq 3U s 


36 R 


— U s ) 

3S 


The circumferential-momentum equation can be stated 


H 3P 

3U G 

3U e U 0 

8U 0 

U 0 Ug 3R 

— •* 

T 0s + T 0r + P H l + 

+ 

Ug 

+ - 

R 36 

3t 

CL? 

CD 

VO 

3S 

R 3S 


Following Hirs' approach [11], the wall shear-stress definitions in 
these equations can be stated 


T S s - — P Us 2 R a mS Cl + (Ue/U s ) 2 ] _2_ 
2 


ms+1 


(5) 


ill 

Tgr - — P Ug 2 R a mr <1 + [(U e -R W )/Ug] 2 | 


nr 

2 

ns 


mr+1 


T0s “ ~ P u s U 0 R a ms Cl + (U e /U s ) 2 ] 
2 


ms+1 


nr mr+1 

T 0r - — p U s (U e - Roo) R a mr {l + C(U e -Ro))/Ug] 2 ( 7 

2 
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where 

R a - 2H U s /v (6) 

The empirical coefficients (ns, ms), (nr, mr) account for different 
surface roughnesses on the stator and rotor, respectively. 

Impeller-Annulus Geometry 

For this analysis, the assumption is made that the impeller is 
nominally centered in its housing. Hence, in the centered position , 
the clearance function H c (s) is only a function of the path coordinate 
S and does not depend on the azimuthal coordinate R6. The inlet- 
clearance function H o (0), the inlet path velocity U s (0), and the inlet 
radius R(0) are denoted, respectively, by Cj_, Vj, and Ri. In terms of 
these variables, leakage volumetric flowrate is defined by 

Q - 2irRiCiVi (7 ) 

The length of the leakage path along the impeller face is defined by 

Zi+L 

L s - / /l+ dR 2 dZ (8) 

Zi dZ 

Nondlmenslonallzatlon and Perturbation Analysis 

The governing equations define the bulk-flow velocity components 
(U s , He) and the pressure P as a function of the coordinates (R6, S) 
and time, t. They are conveniently nondimens ionali zed by introducing 
the following variables 

u s - U s /Vi , U0 *= U0/Rju>, p - P/pVi 2 

h - H/C if s - S/L s , r - R/Ri (9) 

t = u>t , b * V^/R^ oj, T * Lg/V j_ 

The objective of the present analysis is to examine the changes in (u s , 
ue, p) due to changes in the clearance function h(e, s, t) caused by 
small motion of the impeller within its housing. To this end, the 
governing equations are expanded in the perturbation variables 

u s “ Ug 0 + eu s i , h - h 0 + ehi 

u 0 ■ u 0o + eu 61 , p - p 0 + epi (10) 

were e-e/Ci is the eccentricity ratio. The following equations result: 

Zeroth-Order Equations 

(a) Path-Momentum Equation 

dp 0 clug 0 1 dr 

+ u s0 (— J 

ds ds r ds 


(1 la) 


2 o s + a r 

(ueo /b ) + ( ) u so 2 “ 0 

2 
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(b) Circumferential-Momentum Equation 
du 0o ueo dr 

2 + 2 + [o r (u 6o ~r) + o s ue 0 ] - 0 (11b) 

ds r ds 


(c) Continuity Equation 

r h c u 30 - 1 (11c) 

observe that the continuity equation follows directly from Eq. (7) 

First-Order Equations 

(a) Path-Momentum Equation 

3pi 

1 


“ h 1 A 1 S " U01 A 2s - u 3l A33 

3s 

3u s1 ueo 9u s1 

- <uT + ojT + u 

3t r 30 

(b) Circumferential-Momentum Equation 
L 3 1 3pi 

b hi A-|0 - UQi A 2 0 - Up-] A30 

Ri r 36 


- uiT 


(12a) 

3u s1 

so " 
3s 


A ie _ u 9i a 29 

_ u s 1 a 39 


3U01 

u 0o 

3u0i 

3u©l 

+ ojT 

+ u so 


9t 

r 

30 

3s 


(c) Continuity Equation 

3ugi o)T 3u0i 1 dr 1 dh 0 

+ — + u s i ( + ) - 

3s r 36 r ds h 0 ds 


^1 u so dh 0 1 

_ . l u so 

hn^ ds ho 


3hi U0 O 3h-| 3h-j 

+ u)T + ujT ) 

ds r 30 3 t 


(12b) 


(12c) 


Most of the parameters of these equations are defined in Appendix A. 
The quantities o s and o r are defined by 

O3 “ (Lq/Hq) Ag , Op “ (L 5 /H 0 )A r (13) 

where A s and A r are dimensionless stator and rotor friction factors 
defined by 


A s ■ ns R ao m ‘ tl + (u0 O /bUg O )^l 
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Zeroth-Order-Equatlon Solutions 

The zeroth-order Eqs. Til) define the pressure and velocity 
distributions for a centered impeller position. For a known volumetric 
flowrate, the continuity equation completely defines u 80 . The 
continuity equations can be used to solve for u 80 to obtain 


du 80 1 dh 0 1 dr \ 

- - u so l + J 

ds h 0 ds r ds 

substituting into Eq. (11a) yields 


dp 0 1 ^dr ^ U0 O 2 

ds r ds b 2 


^o r +c s ^ 1 dh 0 1 dr 

2 h n ds r ds 


“so 2 


du 0o U0 O dr 

+ + [o r (u0 O -r) + o s U0 O ]/2 - 0 

ds r ds 


<m> 


Those equations are coupled and nonlinear and must be solved 
iteratively. The initial condition for U0 O (O) is obtained from the 
exit flow condition of the impeller. The inlet and discharge pressure 
of the impeller are known and serve, respectively, as the exit (P e ) and 
supply (P 8 ) pressures for the leakage flow along the impeller face. 
The inlet conditions for p Q is obtained from the inlet relationship 

P 8 - P 0 (0, 6, t) - p (1+5) U so 2 (0, 6, t)/2 (15) 

From this relationship, the zeroth-order pressure relationship is 

P o (0) - P s /pVi 2 - (1+0 u so 2 (0)/2 (16) 

The Impeller exit may also include a restriction yielding a 
relationship of the form 


P 

P (L 8 ,6,t) - P e - " C de U s 2 (L s ,0,t) (17) 

2 

The solution to the zeroth-order Eqs. (14) must be developed 
Iteratively since all of the coefficients depend on the local path 
velocity U 80 . In this study, the equations are solved by the following 
iterative steps: 

(a) Guess or estimate which then defines u 80 (s). 

(b) Calculate p o (0) from Eq. (16), and use a specified uq o (0) 
as initial conditions to numerically integrate Eqs. (14) out 
to s - 1 , i.e., the annulus exit. 
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(c) Based on the difference between a calculated exit 
pressure and the prescribed exit pressure, calculate a 
revised and repeat the cycle until convergence is achieved. 


First-order Equations Solutions 
The first-order Eqs. (12) 
u s1 (s, 6 , t) , uei (s,6,t), and pi 
clearance function h-j . From Eqs. 


define the first-order perturbations 
(s,0,t) resulting from the perturbed 
(H) and (9), h-j can be stated 


eh! - { - 


x+av ( — ) z 
Ci 


.Ri 


cosY + ay ( — ] rsinY } cos0 
0 i 


f- 


y-a x (— ) z 
Ci 


cosY - a x ( — ) rsinY J sine 


( 18 ) 


- hi c (s,T)cose + hig (s,t) Sine 

The theta dependency of the dependent variables is eliminated by 
assuming the following, comparable solution format 


U s1 “ u s1c cos 6 + u s1s sine 
uqi - UQi c cose + uqi s sine 
Pi “ Pic cos © + Pis sine 


Substituting into Eqs. (12) and equating like coefficients of cose and 
sine yields six equations in the independent variables s,t. By 
introducing the complex variables 


Usl - u 3lc + J u s1s , U01 - UQIc + J u ei3 

£l “ Pic + JPls * £l “ h 1c + J h 1s 


(19) 


these real equations are reduced to the following three complex 
equations in s and t. 


9£i_ 

3s 


Hi A 1s " £91 a 2s “ £s1 a 3s 


<* u s1 U0o Gugi 

wT — " j wT Ug-] + u P q 

8t r 8s 


b L e 


- J - (“) £i - hi *10 - U0i A 20 _ £s1 a 30 

r Rl 


3u 0 i U0 O 3UQ! 

u)T — j ooT _U 0 1 + UpQ 

3t r 3s 


(20a) 


(20b) 
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9 iis 1 

uT 

1 dr 

1 dh 0 


- J “ U 01 + Usl 

f— — 

+ — 

3s 

r 

r ds 

h 0 ds 


hi 


u so dho uT u 0o 


h 0 2 ds 


h 0 r 


From Eqs. (3) and (17), hi can be stated 
ehi - -q 


L dz 

( — ) — “ <* G o 


L s ds 


u so «|hl <*> T 3hl 

h^ 3 s hg 3 t 


where 

L 2 dz Ri 2 dr 

G 0 - ( ] z — + ( ) r — 

C^Lg ds C^Lg ds 

q - x+Jy , a - ay - Ja x 

From Eq. (21), the following additional result is obtained 

3hi L d 2 z 

e - ~ q ( — ) — ~ - a F-| 

3s L s ds 2 


(20c) 


( 21 ) 


(22) 


(23) 


where 


L 2 

dz 2 d 2 z 

Hi 2 

dr 2 d 2 r 

) 

( — ) * z — 

- ( ) 

(-) +r -_ 

C i Ls 

ds ds 2 

C i L g 

ds ds 2 


The time dependency of Eqs. (20) is eliminated by assuming 
harmonic seal motion of the form 

JfT jf T JfT 

q - q Q e a - a c e hi - hi c e (25) 

f - n/u> 


where fl is the seal whirl frequency and q Q and a 0 are real constants. 
The associated harmonic solution can then be stated 


JfT JfT 

Hsl ” usi e U 01 - Tieie 

JfT 

Pi - Pi e 


(26) 
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Substitution from Eqs. (25) and (26) into the governing complex partial 
differential equations yields the following three complex ordinary 
equations in s 


d 

ds 

where 
[A] = 



and 

r 


f 2 

f 3 



+ [A] 







(27) 

( 28 ) 


1 dr 1 dh 0 

r ds h 0 ds 

a 30 /u so 


wT 

-j - 0 

r 

b Lg 

(A 2 e + jrT)/u so -j (— 

ru 30 Ri 


1 dr 1 dh Q 

a 3s " u so ( + ) + JfT A 2s ♦ juTugo 

r ds h Q ds 


TT dz 
F 2 + J 

I ho ds 

A i g dz 

u 50 

dz TT dz 

“ A is " u so^2 - J u so 

ds h 0 ds 


0 


(29) 


F 3 + jG 0 rT/h 0 
_A 10^O' / ' J SO 

-G 0 A l s ~ u so f 3 - JGq u so n/h 0 


(30) 


= w (f -U0 O /r) 

u so z 1 dh Q dz 

ho ds 2 h 0 ds ds 



ho 


Go 

Ho 



ds 


(3D 



The following three boundary conditions are specified for the solution 
of Eq. (27) : 


(a) The entrance-perturbation, circumferential velocity is zero, 
i.e. , 

ue -1 (0) - 0 (32. a) 

(b) The entrance loss at the seal entrance is defined by Eq. 
(15), and the corresponding perturbation-variable 
relationship is 

pi (0) - - (1 ♦ O usl (0) (32. b) 

(c) The relationship at the exit is provided by Eq. (17) and 
yields the following perturbation relationship 

Pi (D - C de u 80 (1) u s1 (1) (32. c) 

The value for C de depends on the wear-ring seal geometry. Solution 
of Eq. (27) for the boundary conditions of Eqs. (32) is relatively 
straightforward, involving successive solutions for displacement and 
rotation excitations. The complete solution is the sum of the 
homogeneous solution (which depends on the unknown initial conditions 
Pi (0), u 9 i (0)) and the particular solutions which are proportional to 
either q Q or a G . Complete solutions are developed separately for the 
two vectors on the right-hand side of Eq. (27), which satisfy the 
boundary conditions of Eqs. (32). By virtue of the problem’s 
linearity, these two solutions may be added to obtain the complete 
system solution or employed separately to calculate the rotordynamic 
coefficients. The solution to Eq. (27), due to displacement 
perturbation, is obtained by setting a Q » 0, and may be stated 



The soluton due to angular perturbations is obtained by setting q 0 - 0 
and may be stated 



^4c + J ^4s 

*5c + J f 5s 
f 6c + J f 6s 



(3*0 


Reaction Forces and Moments 

From figure 5, the differential force components acting on a 
differential-impeller surface area can be stated 


dFx - - (PcosY - x r «, sinY) R d6dS cos© + x r Q sin© R d©ds 

dF Y - - (PcosY - t r9 sinY) R d6ds sin© - x r g cos© R d©ds (35) 

dFz - - (PsinY + x r9 cosY) R d6ds 
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The shear-stress contribution to these differential-force components 
ha? been neglected in prior analysis of seals. The X and Y components 
are used to define rotordynamic coefficients; the Z component 
defines the axial thrust. Taking moments of the differential force 
vector about the origin of the X, Y, Z system yields the following 
differential moment components. 

dM x = { -[P(RsinY-ZcosY) + T rs (RcosY + ZsinY)] sine + Z T r QCOs6 } 

R d9ds . 

dMy = { -[P(RsinY-ZcosY) + -t rp (RcosY + ZsinY)] cosG + Z x r0 sine } 

R d6ds " (36) 

dM x = “ t r Q R 2 d0ds 

The X and Y components yield rotordynamic coefficients; the Z component 
defines the drag torque. 

From Eqs. .(35) and (36), the force and moment perturbations are 
stated 

L s 2 tt 

F X1 = ~d" / (Pi cosy - x r3 i sinY) cos0 R d6dS 
o o 

(37a) 

L s 2 ir 

+ ef / Tr91 9 i n 9 R d9dS 
o o 


L s 2ir 

Fyi = -el f (PicosY - x rs i sinY) sine R dddS 
o o 

(37b) 

L s 2 it 

-e/ / T r Qi cosG R d6dS 

o o 


L p 2 tt 

M X i = -el" S [Pi (RsinY - ZcosY) + t pp i (RcosY+Z sinY)] R sine dedS 
o o 

(37c) 

L s 2ir 

+e/ / i r 0 i ZR co30 dedS 

o o 


Ijg 2 IT 

Myi = el f [Pt (R sinY - Z cosY)+ T rs i (RcosY+Z sinY)]cose R dedS 
0 o 

(37d) 

Lp 2 -rt 

+ ej" I x r 0 i ZR sine dedS 

o o 
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From Eq. (5), the perturbation shear stresses can be stated 
T sr 1 = ( B s 1 u s i + B S 2 uqi + Bg^ hi ) 

T 0 r 1 = pVi 2 (Bei Ug i + Bq 2 uqi + BQ3 hi ) 


(38) 


The coefficients of these equations are defined in the appendix. 
Successive substitutions from (a) Eqs. (9) and (18), (b) Eqs. (19), and 
(c) Eqs. (21), (24), and (26) into Eqs. (37) yields 


-jft 


F r + J F 0 ( F X1 + J f Y 1 ) e 


(39a) 


eirL< 


CdL 


( L 

dz 

’ R dr 

7 Pi (_) _ + Us1 

(— ) — Bgi + j Bqi 

l Lg 

ds 

L g ds 


R dr 

U 01 |( ) b s2 + j b 02 

Lc, ds 


hi 


R dr 

( ) Bg 3 + j Bq 3 


L s ds 


} 


rds 


M 0 - J M r (Myi-jMxi)e 


-jft 


FoL 


where 


Note that 


F 0 L 


eirLg 1 ( 

,» 1 2 , Or 

L 

dz 1 

— / J*1 

( ) i + 

(- z - 

c d 0 l 

LLg ds 

L s 

ds i 


+ U s 1 


+ U 01 


Ri dr dz 

( — ) (z r ) B 3 i + j Bqi 

L,g 


ds ds 


Rj dr dz 
(— ) (z r ] B s2 + j B 02 


ds ds 


hi 


Ri dr dz 

( — ) (z r ] B s3 + j B 63 


ds ds 


I 


rds 


F 0 = 2 RiLAP 


(39b) 


(40) 


AP = Pg“P e -C(j 


pVi^ 


(41) 
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i? the total pressure drop along the leakage path from Impeller 
discharge to inlet. 


Rot ordvnamlc Coefficients for Displacement Perturbations 

Eqs7T39) apply for simultaneous displacement and slope 

perturbations and are solved alternately for displacement and slope 
perturbations. For the displacement solution, a 0 is set equal to zero, 
and Eqs. (21) and (33) yield 


F r (f) 

.o 

H L s 1 r L dz 

- — (-) / f3c (“) - 

C d L o L Lg ds 

f rq " “ 

Qo f o 


R a dr 

+ f£ C ( — ) B s2 “ f 2s E 

L s ds 

F 0 (f) 

f 0q “ 

Qo f o 

it 1 T L dz 

- - (-) / f 3 » (-) - 

C d L o L L s ds 


R^ dr 

+ f“2s — b s2 + f 2c 1 

Lg ds 

Mq ( f ) 

“Oq “ 

q 0 f o l 

IT Ls 1 r 

“ ( — ) S f 3? E, + f 1c 1 
C d L o L 


L dz 

_ f 2s b 62 " ( ) e 2 

Lg ds 

M r (f ) 

it Ls i r 

— (-) / f 3c E-| + f 1s 

C d L o L 

m rq * ® 

Qo^oL 


L dz 

+ f 2c b 02 " ( — ) B 

Lg ds 

where 

R t 2 dr L 2 d 

Ei - ( ) r — + (— ) z - 

LLg ds Lg d 


R< dr dz 

E 2 . (-) U - - r -) 

Lei ds ds 


LRi dz dr 1 

(— ) b.3 

L 2 s ds ds J 


rds 


(42) 


rds 


rds 


(43) 
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The right-hand side of Eqs. (42) is only a function of the 
frequency-ratio, f, and can be evaluated for selected values of this 
parameter. 


The comparable results from Eq. (1) are 


F r (f) 

f r q(f) “ - - 

Qo f o 

fe q (f) ■ - 

Qo F o 

M 0 (f) 

ni 9q(f) “ - “ 

Qo f o^ 

Mp(f ) 

-m r q(f) - - 

QoFqL 

Where the nondimensional 


(K + fc - f 2 M) 
k - f C 

(K + fc - f 2 M ) 
ae ae ae 

K - f C - f 2 m 
ae ae ae 

coefficients are defined by 


(44)' 


K - KCi/F 0 , k - kCi/F 0 , fl - MCi oo 2 /F 0 
C “ CCj 0)/F o , c ■ cC^to/Fr, 

. 1 . (45) 

K - K Ci/LF 0 , C - C C. ui/LF ot M - M C j( d 2 /LF 0 

ae ae 1 v ae ae i ° ae ae i ° 

k - k Cj/LF 0 , o - c C.a)/LF 0 , m C. oo 2 /LF 0 

ae ae 1 v ae ae i u ae i ° 

The dynamic coefficients are obtained by equating the right-hand sides 
of Eqs. (42) and (44) and by carrying out a least-square curve fit on 
the results from Eq. (42). 


Rotordynamlc Coefficients for Slope Perturbations 

Applying the procedure of the preceding section to Eqs. (42) with 
q 0 “ 0 yields 



F r (f) 

ir 

L 

1 r , L 

dz 

Ri 

dr 




fr 

(f) - 

m — ■ 

(-) 

f Uec (-) 

— + file 

(~) 

— 

B s1 

flJS 

B ei 


«o F o 

Cd 

Ls 

o L Ls 

ds 

Ls 

ds 







Ri 

dr 


Ri 

dr 


1 




+ f 5c 

(-) 

— B s2 - f5s 

®62 " G 0 

(-) 

— 

b s3 

rds 





Ls 

ds 


Ls 

ds 

“ 

1 



F e(f) 

ir 

L 

1 r L 

dz 

Ri 

dr 




f e 

(f) - 

- - — 

(-) 

f I f 6s ( ) 

— + fils 

(~) 

— 

B s1 

+ f 4 c 

B ei 


a o F o 

Cd 

Ls 

o L Ls 

ds 

Ls 
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rds 


Ri 

+ f 5s b s2 + f 5c b 82 

Ls 



M 0 (f) 

m^a(f) - 

ocqFoL 


IT L Ip 

(— ) / f 6c E i + r*c e 2 B si “ f 4s B ei 

C d Ls o L 


+ f5c e 2 b s2 - *5c b 82 “ G 0 E 2 B s3] rds 
M r (f ) It L 1 r 

m r a ( f ) - ( — ) / I f6s E 1 + f 4s e 2 b s 1 + f 4c ^1 

a 0 F 0 L C d Ls o 

+ f 5s E 2 B s2 + f 5c Bq 2 “ G 0 B e3 ] rds 

Eq. (1) yields the corresponding solution format 


( 46 ) 


W f ) 

m 9a ( f ) 


F r (f ) 

- - (Kea + f cea - f 2 Ma) 

a o F o 


F 0 (f) 

a o E o 
Mg ( f ) 

a 0 F 0 L 


■ (k - fC - ) 

ea ea ea 


- (k + fc - f 2 M ) 
a a a 


(47) 


~ni ra ( f ) 


M r (f) 

a 0 F 0 L 


k - fC - f m 
a a a 


where 


K - K /F 0# k - k /F 0 , M - n w/r 0 
ea ea ° ea ea ° ea ea ^ 

C * C oj/F 0> c *= c u>/F 0 m = 1110 ) /F 
ea 0> ect ea u * ea ea 0 


ea 


(48) 


u)^/F f 

K - K /F 0 L, k - k /F 0 L , M « M uj 2 /F 0 L 
a a 0 a a ° a a u 

E a “ C a */F 0 L, c a - C a a)/F 0 L, = m a w 2 /F o L 

The rotordynamic coefficients are obtained by equating the right-hand 
side of Eqs. (46) and (47). 


Predictions and Comparison to Experimental Results 

Figures 6(a) and 6(b) illustrate, respectively, a nominal pump- 
impeller geometry with a conventional wear-ring seal and a modified 
pump Impeller with a face seal. Bolleter et al. tested the face-seal 
impeller to eliminate the forces which would normally be generated by 
the wear-ring seal. Their tests were at best efficiency point (BEP) 
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with the pump running at 2000 rpm, while developing 68m of head and 1 30 
1/sec of flow rate. The impeller has seven blades and an impeller exit 
angle of 22.5°. The test fluid is water at 80°F. 

The present analysis requires an estimate of the AP across the 
impeller versus the total head rise of the stage. At U. Bolleter's 
suggestion, the impeller AP was estimated to be 70? of the total AP of 
the pump. An estimate of the inlet tangential velocity is also 
required. Fortunately, pitot-tube measurements are available, 
indicating that the inlet tangential velocity is approximately 50? of 
the exit impeller surface velocity; i.e., uq 0 (o) ~ 0.5. This is in 
contrast to a theoretically-predicted tangential velocity of 0.72 RjO> 
based on the blade-exit angle. Adkins' measurements at Cal Tech [12] 
also show substantially lower exit tangential velocities than predicted 
from exit blade angles. 

Both walls of the impeller were assumed to be smooth and 
represented by Yamada’s [ 13 ] test data; mr - ms - -0.25, nr - ns - 
0.079. The inlet loss for the impeller, £, was assumed to be 0.1. The 
discharge coefficient for the seal was calculated Iteratively as 
follows. With an assumed C de , Eqs. (HO, (15), and (17) were used to 
calculate the leakage through the impeller annulus and the pressure and 
tangential-velocity upstream of the seal. The seal is then analyzed 
(with the same equations) using the calculated seal inlet 
pressure and tangential velocity to determine leakage and C de . The 
iteration continues until the leakage predictions for the exit seal and 
the impeller annulus agree. 

Figures 7(a) and (b) illustrate the predicted radial and 
tangential force coefficients f r q and f$q versus the whirl frequency 
ratio f = q/u) for the face-seal, impeller. Results are presented for 
U 0 O (O) “ 0.5, 0.6, and 0.7. The uq o (0) «= 0.5 data of these figures is 
generally consistent with expectations based on experience with seals 
except for a slight "dip" in f r q and "bump" in fgq. However, the peaks 
exhibited at higher value for u@ 0 are quite unexpected. They arise 
primarily due to the centrifugal acceleration term in the path-momentum 
equation. If the term 

2u eo dr 

/ b 2 , 

R ds 

is dropped from the k 2 s definition of Appendix A, the "peaks" are 
substantially eliminated from the force predictions. 

Figure 8(a) and (b) provide predictions for f rq and f 6 q for the 
conventional wear-ring-seal impeller of figure 6(a). The predicitons 
are only for the impeller and do not include the exit wear-ring seal. 
The f Q q results for the two impellers are quite similar; however, the 
f r q values are generally larger for the face-seal impeller. 

Table 1 below provides zeroth-order-solution results for the 
conventional-seal impeller of figure 6(b). 
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Observe that the leakage is reduced by increasing uq o (0). Also observe 
the relatively high seal-inlet-tangential velocity prediction for the 
exit seals, which will predictably lead to increased cross-coupled 
stiffness coefficients and decreased rotor stability. The Reynolds 
number values suggest that care should be taken in extrapolating 
impeller force data from one operating condition to another. 
Specifically, changes in temperature or operating media would be 
expected to yield significant changes in force data. 

The freqency-dependency of f rq and fg q exhibited in figures 7 and 
8 for uq o (0) = 0 . 6 , 0.7 can not be modeled by the rotordynamic- 
coefficient model on Eq. (1). Stated differently, The quadratic 
dependency of f rq and f eq on f, which is specified in Eq. (4M), is 
simply not true. A significantly more complicated dependency is 
clearly in order. While these results of figures 7 and 8 are 
surprising, recent data from the Cal Tech program by Franz and Arndt 
[1-411 are qualitatively similar. 

The U 0 O (O) = 0.5 results of figure 7 are reasonably modeled by a 
quadratic dependency of f and can be modeled by rotordynamic 
coefficients. A comparison of predicted and measured coefficients is 
provided in Table 2. 

Keeping in mind that the present theory does not account for the 
momentum flux exiting from the impeller or the pressure forces 
on the impeller exit, the comparison between theory and experiment of 
Table 1 is encouraging. The prediction of C and M are good. The 
results for k are consistent with Adkins [12] statement that the 
impeller annulus accounts for approximately one half of the measured 
stiffness values in Cal Tech test results. The results for K and c are 
obviously dissappointing; perhaps these coefficients depend more 
heavily on the impeller-diffuser interaction forces. 

Conclusions and Extensions 

An analysis has been developed for the forces on the shroud of an 
impeller. The "bulk-flow" nature of the analysis restricts its 
applicability to impellers having fairly small clearances between the 
impeller shroud and casing. 

The results of the analysis are unexpected in that resonances of 
the fluid system are predicted at inlet tangential velocities which are 
higher than approximately 0.5 Riw. Conventional rotordynamic- 
coefficient models are not an adequate representation of rotor forces 
if the resonances are present. The resonance phenomenon predicted by 
this analysis are in qualitative agreement with recent measurements 
from Cal Tech [14]. 

Rotordynamic-coeff icients predictions from this analysis are in 
reasonable agreement with test results from Bolleter et al. [7], for 
the direct damping and cross-coupled stiffness coefficients, C and k. 

The analysis which was developed in this paper for impeller 
shrouds can also be applied to seals and provides an expanded 
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capability with respect to clearances as a function of axial or path 
coordinate. Prior analyses have been developed for either constant- 
clearance or convergent-tapered seals. The present analysis would 
apply for any continuously varying clearance function. This option can 
be directly used to account for elastic deformation of the seal bore 
due to pressure or other loading. It also provides possibilities for 
optimizing the clearance function to maximize stiffness, damping, etc. 

The addition of shear-stress contributions in the present 
analysis made very little difference in seal calculations. As yet, no 
comparisons have been made on impeller force calculations with and 
without the shear-stress contributions. 

The importance of the moment coefficients, which are defined by 
the represent analysis, are as yet undetermined. Current rotordynamic 
codes are not yet up to accept these coefficients; however, a code is 
under development at TAMU which will account for these coefficients and 
can be used to determine their importance. 

The significance of the "resonance" phenomenon, which is predicted 
by the present analysis, on rotordynamics also remains unsettled. An 
analysis of the phenomenon is planned which will include curvefitting 
the f r q(f) and fgq(f) predictions and examining their influence on a 
Jeffcott rotor model. 

The present analysis can be readily extended to account for axial 
force coefficients. An extension of this nature is projected for the 
coming year. 


APPENDIX A 


Perturbation Coefficients 

A 1s “ [o s (1-ms) + o r (1-mr)] u so 2 /2h 0 
2u 6o dr _ 

A 2s /b 2 + [o r (mr+1) B 0 + o s (ms+1 ) ] u PO /2 

r ds 

d u S o 


a 3s 


2A ie 


2A 20 


+ [(2 + mr) o r + (2 + ms) o s ] u 30 /2 

ds 

-[ ( 1 +mr) Op B 0 (ue 0 -r) + (1+ms) o 3 Bi U0 O ] 


u 30 [(1-mr) ( UQ 0 -r ) o r + (1-ms) u$ Q o«, + 2 


u so (op+o p ) + Op (mr+1) (u 6o -r) 6 0 + 2 

r 

+ o 3 (ms+1) U0 O Bi 


/ 2 

^0o 

r 

dr 

ds 


dr 

— ]/ho 
ds 
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2A^e « o r (u 0 O -r) [mr - (1+mr) Bo (u0 O -r)/u 3O ] 
+ °s u 6o tms - (1+ms) B 1 u^ 0 /u ?0 ] 

B 0 - (u eo -r)/b 2 Uc , 0 {1 + [(ue 0 -r) /bu ?0 ] 2 } 

01 “ ue 0 /b 2 u ?0 [1 ♦ (u 0o /bu so ) 2 ] 


T r0 Perturbation Coefficients 


Bqi - A r (1+mr) (u0 O -r)[1-B o ( u eo _r ) /u so^ /2b 

b 02 - * r [u 50 + (1+mr)(u0 O -r)Bo]/2b 
B 0 2 = Ip mr (u 0o -r )u 30 /2bh 0 
Tp a Perturbation Coefficients 


B s1 ■ ^ r C( 2+mr ) u so - O+mr) Bo ( u 9o _r )3 /2 
Bg2 ” Xp(1 + mr) 8 q ^so^ 2 
b s 3 “ ^r m r u so^^ 2b o 
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U eo (0)/Ri a) 

0.5 

— 

0.6 

0.7 

U 0O (1 )/R(1)« 

0.884 

0.947 

0.996 

m(Kg/sec) 

4.03 

3.75 

3-44 

Rao * 2ViCi/v 

73,300 

__ —— 

68,100 

62,600 


Table 1 . Zeroth-order- solution results for the 
conventional -wear-ring impeller. 



Measured 

Theory 

Face-Seal 

Impeller 

K(N/m) 

-.5*106 

-.042*10 6 

k(N/m) 

.6x106 

.288x1 0 6 

C(N sec/m) 

2570* 

2020 

c(N sec/m) 

7610 

2290 

H (kg) 

29.6 

8.96 

m (kg) 

10.8 

-0.009 


♦Combined viscous and hysteretic damping. 

Table 2. Theory versus experiment for the face-seal impeller. 
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IMPELLER SHROUD 



Figure 1. lapel ler stage. 



Figvire 2. Impeller surface geometry. 



Figure 3- Local attitude angle of impeller surface 



Figure 4. Differential-fluid element. 




Figure 6. Nominal configuration of Sulzer teat Impeller [7] 
with conventional wear-ring seal. 
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FACE- SEAL IMPELLER 



-1.25 -.75 -.25 .25 .75 1.25 1.75 2.25 





-1.25 -.75 -.25 .25 .75 1.25 1.75 2.25 




Figure 7. Nondlmensional force coefficients for the face-seal 
impeller; A tangential-force coefficient, 

B radial-force coefficient. 
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1.25 -.75 -.25 . 25 . 75 1.25 1.75 2.25 

f e n /« 



-1.25 -.75 -.25 . 25 . 75 1.25 1.75 2.25 

f = 


Figure 8. Nondlmensional force coefficients for the conventional 
impeller; A tangential-force coefficient, 

B radial-force coefficient. 
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